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Abstract 



The asymmetries between charm and anticharm mesons observed in fixed- 
target photoproduction experiments are an order of magnitude larger than 
the asymmetries predicted by conventional perturbative QCD. We show that 
these charm meson asymmetries can be explained by a heavy-quark recom- 
bination mechanism for heavy meson production. In this process, a charm 
quark combines with a light antiquark from the hard-scattering process and 
they subsequently hadronize into a state including the charm meson. This 
recombination mechanism can be calculated within perturbative QCD up to 
some nonperturbative constants. After using symmetries of QCD to reduce 
the number of free parameters to two, we obtain a good fit to all the data on 
the asymmetries for charmed mesons from the E687 and E691 experiments. 
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I. CHARM-ANTICHARM ASYMMETRIES 



High energy charm photoproduction is an excellent testing ground for our understand- 
ing of strong interactions The charm quark mass provides a large scale justifying the 
use of perturbative QCD. Theoretical uncertainties in photoproduction are smaller than in 
hadroproduction because only one hadron is present in the initial state. Charm events also 
constitute a higher fraction of the total hadronic cross section than in hadroproduction, 
which facilitates experimental studies. 

High energy charm photoproduction experiments 0-f| have observed an interesting pro- 
duction asymmetry between charm and anticharm mesons. The E691 |2j and E687 || 
experiments at Fermilab measured the cross sections for various D and D mesons from a 
photon beam on a beryllium target. The experiments imposed a cut on the Feynman xp 
variable, xp > 0, which corresponds to the forward direction of the photon beam in the 
7 A/" center-of-mass frame. E691 and E687 found statistically significant asymmetries in 
D+-D-, D°- D and D* + - D*~. In each case the cross section for the anticharm meson 
is greater, though the size of the asymmetry is different for each species of meson. Produc- 
tion asymmetries for and A+ are also measured, but with much larger errors so they 
are consistent with zero. For D + and D~ mesons, the E687 experiment has also measured 
the dependence of the charm asymmetry on the photon energy, E 7 , the meson transverse 
momentum, p±, and Feynman xp. 

The observed charm asymmetries in photoproduction are an order of magnitude larger 
than the asymmetries between charm and anticharm quarks predicted by perturbative QCD. 
Thus they cannot be understood within the fragmentation mechanism for heavy meson pro- 
duction which is based on the factorization theorem for inclusive single particle production 
||. In this mechanism, the c and c are produced in a hard-scattering process and then 
fragment independently into charmed hadrons. At leading order in a s , charm is produced 
by photon-gluon fusion |J which produces c and c symmetrically. The next-to-leading order 
(NLO) correction |7]-|9| increases the normalization of the total cross section by about 50% 
at the energies relevant to the E687 and E691 experiments. The shapes of the xp and p± 
distributions are similar to those at leading order, but there is a tiny charm asymmetry in 
these kinematic distributions, arising from the processes "f+q(q) — * c+c+q(q). For E 1 = 200 
GeV, the prediction from a NLO calculation for the ratio of the c and c cross sections in the 



forward region is about R = 1.006 Jl0|. The asymmetry, (1 — R)/(l + R) = —0.003, is about 
an order of magnitude smaller than the asymmetries for D mesons measured by the E687 
and E691 experiments. Furthermore, the NLO calculations cannot explain the differences 
in the asymmetries of D mesons with different light quark flavors. It is commonly assumed 
that a nonperturbative hadronization mechanism involving the remnant of the nucleon or 
photon after the hard scattering is responsible for the charm asymmetries. 

In this paper, we argue that the charm asymmetries can be explained naturally and 
economically by the heavy-quark recombination mechanism for heavy meson production in- 
troduced in [PH| . In the cq recombination mechanism, a light q from the target nucleon 
participates in the hard-scattering process that produces the c and c. The light q emerges 
from the hard-scattering with small momentum in the c rest frame, and the c and q subse- 
quently hadronize into a final state including the D meson. There is also a cq recombination 
process in which a light q from the target nucleon recombines with a c to hadronize into a 



D meson. Because the parton distribution functions of u and d in the nucleon are greater 
than the parton distribution functions for u and d, the recombination cross sections for cu 
and cd are greater than for cu and cd. This gives rise to the asymmetries between charm 
and anticharm mesons. Since the photon participates in the hard scattering and can couple 
to the light quark, the recombination cross sections for processes involving u quarks differ 
from those involving d and s quarks because of their different electric charges. Thus, the 
heavy-quark recombination mechanism can account for the differences in the asymmetries 
of D mesons with different light quark flavor quantum numbers. Below we will see that it 
can also adequately describe the kinematic dependence of the asymmetries on E 7 , and 

Xp . 

Before discussing the heavy-quark recombination mechanism in detail, we briefly review 
the models which have been proposed in the literature to explain the D meson asymmetries. 
Though these models differ in details, all generate the asymmetries by a nonperturbative 
hadronization process that involves the remnant of the nucleon or photon after the hard 
scattering. The predictions of these models depend on unknown and ad-hoc functions that 
specify the distribution of partons in the remnant. 

The model for the asymmetry that has been most thoroughly developed involves the 
fragmentation of strings connecting the c and c to the nucleon remnant [|12j. After the c 
and c are produced via photon-gluon fusion, the remnant of the target nucleon is left in a 
color-octet state. The remnant is modelled as a color-triplet "bachelor quark" and a color- 
antitriplet diquark with momentum fractions x and 1—x specified by a distribution function, 
F(x). The subsequent hadronization is calculated using the Lund string fragmentation 



model Color-field strings are introduced between the diquark and the c and between 
the bachelor quark and the c. The prediction for the asymmetry depends sensitively on the 
choice of F(x) as well as the relative probability for each string to break into a charmed 
meson or a charmed baryon. This model reproduces the qualitative shapes of the kinematic 
distributions for the D + — D~ asymmetry measured by the E687 collaboration M. However, 
the magnitude of the asymmetries is sensitive to the unknown function F(x). For a "soft" 
distribution that diverges as 1/x at small x (the default option in the PYTHIA 5.6 event 
generator program), the model tends to overestimate the asymmetry in all kinematic regions. 
However, the data can be fit with a "hard" momentum distribution proportional to 1 — x [§J. 

Models for the asymmetry that involve the remnant of the hadronic component of the 
photon have also been proposed. In one such model fp70| , a c is produced by the hard 
scattering process qq —>■ cc involving a q from the photon. The c then recombines with a q 
from the photon remnant to form a D meson. The predictions depend on the distribution 
function for the q in the photon remnant and on a recombination function that could be 
calculated in principle using the Lund string fragmentation model. Another model is based 
on the possibility of intrinsic charm |14| in the photon. In this model ||15|| , a q in the resolved 
photon undergoes a hard scattering and the D forms from a c and q in the photon remnant. 
The prediction depends on the distribution function for the momentum fractions of the c and 
q in the photon remnant and a recombination function that depends on several variables. 
Both of the models involve arbitrary functions that can be tuned to fit the data. 

In contrast with these models, the heavy-quark recombination mechanism |Tl| generates 
the asymmetries within a hard-scattering process. The asymmetries are calculable in per- 
turbative QCD up to some nonperturbative parameters related to the probability for the c 
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and q produced in the short- distance process to hadronize into a D meson. Because the c 
and q do not necessarily have the same color, angular momentum or light quark flavor quan- 
tum numbers as the final state D meson, several nonperturbative parameters appear in the 
calculation. However, it is possible to use heavy quark symmetry, SU(3) flavor symmetry 
and large N c arguments to reduce the number of parameters to two. The asymmetries can 
be adequately described by fitting these two parameters. Because of this, the heavy-quark 
recombination mechanism gives a very economical description of the observed asymmetries 
and does not rely on model-dependent assumptions about nonperturbative hadronization. 



II. HEAVY-QUARK RECOMBINATION MECHANISM 

The factorization theorems of QCD |J show that the leading contribution to the D 
meson photoproduction cross section at large transverse momentum can be written as: 

da[y + N + X]=Y,fi/N® da[y + i -> c + c + X] ® D-_^ , (1) 

i 

where fam is the parton distribution function for the parton j in a nucleon and D-_^ is 
the nonperturbative fragmentation function for c to hadronize into D. At leading order 
in a s , the only parton that contributes is the gluon, so photon-gluon fusion dominates the 
cross section. At large p±, Eq. ([]]) receives corrections suppressed by powers of Aq C d/p± or 
m c/p±- The heavy-quark recombination mechanism introduced in [|ll|] is an 0(AqcD m c/ Pj_) 



power correction to Eq. ([!]). The cq recombination contribution to D production can be 
expressed as 

rfa[ 7 + N - D + X] = f q/N <g> d «h + Q ~> + c] p[(cq) n -> D] . (2) 

n 

Making the replacement c «-> c, q — > q and D — > D in Eq. @ gives the cq recombination 
contribution to D production. The notation (cq) n in Eq. (|2|) indicates that the q in the 
final state has momentum of O(Aqcd) in the c rest frame, and that the c and the q are in 
a state with definite color and angular momentum quantum numbers specified by n. The 
factor p[(cq) n — > D] is a nonperturbative constant proportional to the probability for (cq) n 
to evolve into a final state that includes the D. The recombination process in Eq. (0) is an 
O^Aqcd^c/PjJ power correction to Eq. flU), because the parton cross section in Eq. (||) is 
suppressed by 0{m 2 c /p\) relative to the parton cross section in Eq. (JTJ) and p[(cq) n — > D] is 
expected to be O^Aqcd/i^c) from heavy quark effective theory scaling arguments. Because 
the momentum of the light quark is O(Aqcd), the production of (cq) n in states with L ^ 
is suppressed by powers of Aq C d/p± or AgcD/ m c relative to S-waves. Therefore, it is only 
necessary to compute production of (cq) n with 1 5'o or 3 Si angular momentum quantum 
numbers. The (cq) n can be in either a color-singlet or a color-octet state. 



Berezhnoy, Kiselev, and Likhoded |16 Hi9|l have developed a model for heavy meson 



production that incorporates some aspects of our heavy-quark recombination mechanism. 
Their model is based upon an extrapolation of the NRQCD factorization formalism for heavy 



quarkonium production j2D|. A D or D* meson is modelled as a nonrelativistic quarkonium 
composed of a c and q with masses m c and m q . They take m q = 300 MeV, consistent with an 
interpretation as a constituent quark mass. Their expression for the cross section is obtained 
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FIG. 1. Feynman diagrams for the recombination process 7 + q — > (cq) n + c. Single lines 
represent a light quark and double lines represents a charm quark. The shaded blob represents 
(cq) n . 



by keeping only specific terms in the NRQCD factorization formula. These terms involve 
S-wave matrix elements, which are distinguished by the fact that their contributions at large 
p± are suppressed only by m q m c /p\. Of the four independent S-wave matrix elements, they 
keep only two: the Color-singlet Model matrix element (Oii 3 Si)) D = 3((9i( 1 S'o)) i:) and the 
color-octet matrix element (Og^Si)) * = 3(08( 1 S'o))' D - Thus their model involves only 3 
parameters: m q , ((9i( 1 S'o)) D and (Os^So)) . In their model, the terms in the cross sections 
that diverge like 1/m 



have the same form as our heavy-quark recombination 
cross sections if we make the following identification between the nonperturbative factors: 



q aS TTlq 



p[cq 1 ( 1 S )^D] 
p[cqsCSo) - D] 



1 



Am 2 n m c 

q 



4m^m r 

q 



(Oi( 3 S 1 )) 



D 



(3) 



Our approach provides a model-independent framework for calculating those contributions. 

19] , contributions corresponding to the NRQCD matrix elements 



In the model of Ref. 



(Cq( 3 Si)) and (O s ( 3 Si)) , which require a spin-flip between the partons and the meson, 
are omitted. In the heavy quarkonium case, these matrix elements are suppressed by v 7 and 
v 3 with respect to (Oi^Sq)) , while (0 8 ( 1 So)) D is suppressed by v 4 . But when the antiquark 
q becomes light, v approaches 1 and there is no parametric suppression that can be used to 
justify omitting them. Furthermore, the NRQCD operators Ci( 1 S'o), Ci( 3 5i), Os^Sq), and 
Os{ 3 Si) mix under renormalization. This implies that at some order in perturbation theory, 
inconsistencies will arise if any of these operators are omitted. It is possible that some of 
them have matrix elements that are numerically small, but they must all contribute at some 
level. 

The most important difference between our approach and the model of ref. [IS] is that 
we use our nonperturbative probability factors p only for very specific terms in the cross 
sections: those that diverge as l/x q , whre x q is the light-cone momentum fraction of the q 



relative to the heavy quark. In contrast, the model of Ref. |19| is applied in all kinematical 
regions. Their parameters m q , (Oi( 1 S , )) D , and (Og^So)) are used to calculate not only 
the heavy-quark recombination contribution to the cross section, but all other contributions 
as well, such as the fragmentation function D c _+d(z). In contrast, we take the fragmentation 
function to be an nonperturbative function that is independent of the probability factors 
p for heavy-quark recombination. We use fragmentation functions which have been been 
extracted from D meson production in e + e~ collisions. With only 3 parameters, the model 



of Ref. JT9[ is very constrained, but the predictive power comes at a sacrifice in rigor. 



The four Feynman diagrams that contribute to the recombination process 7+g — > (cq) + c 
at lowest order in a s are shown in Fig. [I]. Since the momentum of the light quark in the final 
state, p' is O(Aqcd) in the c rest frame, its transverse momentum can be neglected and 
its light-cone momentum fraction, x q , is 0(AQCD/ m c)- in the limit x q — > 0, the diagrams in 
Figs. 0a), b) and d) are 0(x~ 1 ) while the diagram in Fig. |l|c) is 0(x q ). Therefore the diagram 
in Fig. |l|c) gives an 0(A<2CT>/m c ) suppressed contribution to the heavy-quark recombination 
rate and can be neglected. The calculation is very similar to the calculation of the process 



g + q —> (bq) + b in Ref. ||11|| . Interested readers can refer to [JTTJ] for a detailed description of 
the calculation. The color-singlet 1 5'o contribution to the cross section is obtained by first 
making the following substitution in the parton amplitude: 

Vi(Pc)uj(p' q ) — ► x q j^-m c f + (0c - m c )7 5 . (4) 

We then set p' = x q Pc in the rest of the amplitude and take the limit x q — > 0. The 
nonperturbative parameter / + is related to the amplitude for the (cq) produced in the short- 
distance process to bind into a D meson. / + is proportional to the D meson decay constant, 
/d, times a moment of the D meson light-cone wave function, and therefore is expected to 
scale as m" 1 / 2 in the heavy quark limit. However, / + does not account for the possibility 
that the (cq) can hadronize into states that include other light hadrons whose momenta are 
soft in the D meson rest frame. In order to take these final states into account, /+ is replaced 
with the parameter pi[cg( 1 S , ) — > D] after squaring the matrix element. Pi[cg( 1 S , ) — > D] is 
proportional to the inclusive probability for the color-singlet, spin-singlet (cq) to evolve into 
any state containing the D hadron. The calculation of the color-singlet 3 Si cross section is 
identical except that 75 in Eq. (f|) is replaced by /, where e M is a polarization vector, and 
f\ is replaced by pi[cg( 3 5'i) — > D]. 

The D production cross section also receives contributions from processes in which the 
(cq) is produced in a color-octet configuration in the short-distance process. The color-octet 
1 5'o cross section is obtained by making the substitution 

Vi(pc)uj(p' q ) — > x Q<^ T ii m cf+(^ - m c )7 5 (5) 

in the parton amplitude, replacing p' — > x q Pc and then taking the limit x q — ► 0. After 
squaring the matrix element and summing over colors in the final state, (/+) 2 is replaced 
with ps\cq( l So) — > D]. The normalization in Eq. (|5|) is chosen so that ps/pi is the ratio of 
the probabilities for a color-octet (cq) and a color-singlet (cq) to hadronize into a D hadron. 
Finally the color-octet 3 Si contribution is obtained by replacing 75 in Eq. ((^) by and 
replacing (/|) 2 by p 8 [cq( 3 Si) -> 25]. 
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The final results for the parton cross sections for color-singlet cq recombination in 7 + q 
collisions are 



fl^'si 1 ')] 



256 ix 2 el a a 2 ml 



+ 



81 S 3 
m 2 c S (_S^_ 2(l + re)£ 



S / kT- 2 

u{ 1 + ^ 



u 2 



T 



k 2 T\ 2m 4 c S 3 
+ 4k + — + 



S 



T 3 U 2 



1 + 



reT 



(6) 
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(7) 



where re = e g /e c is the ratio of the electric charge of the light quark q to that of the charm 
quark, S = (p g + p 7 ) 2 , T — (p 7 — p c ) 2 — m 2 and U = (p 7 — pc) 2 — wic- These variables are 
defined so that S + T + f/ = 0. If ^ is defined as the angle between the incoming q and the 
outgoing (cq) in the parton center-of-mass frame, then T = -S/2[l - (1 - 4m 2 J S) 1 / 2 cos 9]. 
The contributions to the cross section for D meson production are obtained by multiplying 
the parton cross sections in Eqs. flfj) and (0) by pi[cq(}So) — ► S] and pi[cg( 3 5'i) — > D], 
respectively. Note that the recombination process is calculated in the heavy quark limit, so 
we do not distinguish between m c and the D meson mass. 

For the photoproduction diagrams in Fig. 0, the color-octet amplitude is identical to 
the color-singlet amplitude up to an overall color factor. The color-octet cross sections can 
be obtained from the color-singlet cross sections in Eqs. (§) and (|7|) by replacing pi with 
P8/8. Therefore, heavy-quark recombination contributions to photoproduction cross sections 
depend on the following linear combination of color-singlet and color-octet parameters: 



Peff = Pi + o PS ■ 



(8) 



It is interesting to compare the relative sizes of the heavy-quark recombination and 
photon-gluon fusion cross sections in different regions of phase space. Setting 9 = it/ 2 and 
expanding in powers of m 2 /S, one finds that the ratio of the parton cross sections is 



da[y + q — > (cq) + c] 



do- [7 + g — > c + c] 



0=71-/2 



64 7r(2-re) 2 m 2 
81 



64vr(2 - re) 
27 



ml 
S 



for cq^So) , 
for cg( 3 Si). 



(9) 
(10) 



For 9 = 7r/2, S = 4:(p 2 ± + m 2 ), so the parton cross section for heavy-quark recombination 
is suppressed by m 2 /p\ at large p±. The recombination cross sections are therefore sup- 
pressed by AqcdttT'c/Pi, relative to the photon-gluon fusion cross section, in accord with the 
factorization theorem. The differential cross section is finite as 9 — > and 9 — > 71 because 
the heavy quark mass acts as an infrared cutoff. Setting 9 = and expanding in rn^/S, one 
finds that the ratio of the two parton cross sections is 
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da[y + q — > (cq) + c] 



dopy + g — > c + c] 



for cqC So) M 3 Si). (11) 



Thus, the heavy-quark recombination cross sections have no kinematic suppression factor 
when the (cq) emerges from the parton collision in the same direction as the incident q. On 
the other hand, the recombination cross sections are suppressed by factors of m 2 c /S when 
the (cq) emerges in the same direction as the incident photon: 



do- [7 + q — > (cq) + c] 



dopy + (7 — > c + c] 



for c^), (12) 

512tt(1 + k 2 ) m\ , 3c , , 10 , 
81 if °^ ^ ' ( ^ 



The kinematic suppression is especially strong in the 1 Sq case. 

To analyze the photoproduction data of E687 and E691, it is necessary to compute the 
cross section for various species of D mesons. It is important to realize that the color, 
angular momentum and light quark flavor quantum numbers of the (cq) n are not necessarily 
the same as the D meson in the final state [p}| . The transition from the (cq) n to the final 
state that includes the D meson is governed by nonperturbative physics that can change 
these quantum numbers. For instance it is possible for the (cq) to be produced in a color- 
octet state in the short-distance process and emit soft gluons during the nonperturbative 
transition to emerge as a color-singlet D meson. It is also possible for the spin of the light 
quark to be flipped during this transition, resulting in a D meson with different angular 
momentum. Finally, a light q'q' pair can be created in the transition, followed by the c 
binding with the q' to form a D meson with different flavor quantum numbers than the 
original (cq) n . Thus for a specific D meson, several nonperturbative parameters enter the 
calculation. For instance, for the D° meson, the parameters are p e fi[cu( l So) — > D°] and 
p e s[cu( 3 Si) — > D°], together with those obtained by replacing u with d or s. 

Symmetries of the strong interaction allow us to greatly reduce the number of free pa- 
rameters. Heavy quark spin symmetry [^TJ implies that the spin of the c quark should 
remain unchanged throughout the nonperturbative transtion. This means that transitions 
which can be related by flipping the spin of the heavy quark in both the initial and final 
state should have the same rate. For example, 

PcslcdCSo) -> D°] = p cff [cd( 3 5 1 ) -> D*°] . (14) 

One can also use SU(3) light quark flavor symmetry to relate nonperturbative transitions 
with different light quark flavors. An example of an SU(3) relation is 

PcslcuCSo) -> D°] = pcsicdCSo) -> D + ] = p eff [cs( 1 5 ) -> D+] . (15) 

Finally, one can argue that in the large iV c limit of QCD, light quark-antiquark pair produc- 
tion is suppressed. This suppression in low energy QCD is the basis for the phenomenolog- 
ically successful Zweig's rule. Suppression of light quark-antiquark pair production in the 
long-distance transition then forces the (cq) n to have the same flavor quantum numbers as 
the final state D meson. 
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Using the symmetries discussed above together with the large N c argument, one is left 
with two parameters. One is 

Psm = Pce[cuCS ) -> D°] = p cS [cdCSo) -> D+] = p eS [csCS ) -> D+] 

= Pos[cuCS 1 ) - D*°] = PcS [cd( 3 S 1 ) -> 7J*+] = pefffc^ 3 ^) - 7J S * + ] . (16) 

The subscript sm stands for spin-matched, since in all these transitions, the spin of the cq 
is the same as the spin of the D meson. The other parameter describes transitions in which 
the light quark has its spin flipped: 

Psf = Pcff [cu{ z S 1 ) -> D°] = p cS [cdCS,) -> D + ] = p^[cs{ 3 Sx) -> £>+] 

= peff[cu( 1 ^ ) - 7T°] = Pcff [cdCSo) -> 7J*+] = Pcff [c^ 1 ^) - 7J S * + ] . (17) 

While the use of heavy quark and 577(3) symmetry greatly reduces the number of parame- 
ters, it should be kept in mind that these relations are only expected to hold at the 30% level. 
It is also important to note that the flavor-changing transitions are not strongly suppressed 
because N c = 3. The relevant nonperturbative parameters for these transitions are expected 
to be smaller than p sm and p s f but nonvanishing. These are not included in the calculations 
of this paper because an adequate fit to E687 and E691 data can be obtained using only 
the two parameters p sm and p s f. However, it may be necessary to include flavor-changing 
transitions in other applications, or if the accuracy of current data on charm asymmetries 
significantly improves. 

It is important to note that when a D is produced by a cq recombination process, 
the associated c quark can produce a D meson via the usual fragmentation mechanism. 
Therefore, the recombination process can contribute to the direct cross section for D mesons 
in two ways: 

a) £ daYi + q^ (cq) n + c] p n [(cq) n - D] , (18) 

n 

b) J2 + 1^ (cq) n + c}p n ® D C ^ D . (19) 

q, n 

In a), the (cq) n recombines into the D meson. In b), the (cq) n recombines into a D meson, 
and the D comes from the fragmentation of the recoiling c. For direct D production, Eq. ([18]) 
and Eq. ( |T9"D are replaced with their charge conjugates. Process a) produces more D~ and 
D° mesons than D + and D°, simply because the nucleon contains more u and d than u 
and d. In process b), JJ's are produced when there is a q in the initial state, while JJ's are 
produced when there is a q in the initial state. The excess of u and d over u and d in the 
nucleon then leads to more D than D. It turns out that in the case of D + — D~ and D° — D 
mesons, process b) dilutes the asymmetry produced in a). For D s mesons, process a) cannot 
generate an asymmetry because the s and the s content of the nucleon are identical. Process 
b) generates an asymmetry for — Dj which has the opposite sign as the D + — D~ and 
D° — D asymmetries. 

Finally, when computing the inclusive cross section for D and D* mesons, one must 
consider the possibility of feeddown from other D mesons that are produced via photon-gluon 
fusion or heavy-quark recombination and then decay to D or D*. The feeddown from excited 
D mesons, e.g. is not included but is expected to be small. Therefore, the inclusive 
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D* cross section is the sum of the photon-gluon fusion and heavy-quark recombination 
contributions, which will be referred to as the direct cross section. The inclusive D meson 
cross section also includes feeddown from D* mesons. Denoting the direct cross sections by 
o"dir an d the inclusive cross sections by <7 inc , one finds 

0.323 <7dir [£>*+] , (20) 
0.677a dir [Zr + ] + cx dir [Zr ], (21) 
^AD*s + ]- (22) 





= cr dlr [ J D + 


cr inc [D } 


= cr dil [D°] 






The coefficients of <7di r [-D* + 


] in Eq. (]2( 



and Eq. (|21|) are the branching fractions of D* + into 
D + and D° ||22|| . The cross sections for the charge conjugate states have similar expressions. 



III. COMPARISON WITH PHOTOPRODUCTION DATA 

In the calculations presented in this paper, the contribution to charm photoproduction 
from the hadronic component of the photon is neglected because it constitutes less than 
5% of the cross section for the fixed-target experiments in the energy range 50 GeV < 
Ej < 400 GeV 0. The charm quark mass is taken to be 1.5 GeV and both the factoriza- 



tion and renormalization scales are set equal to J ml + p\. We checked that varying this 



scale between Jm^ + p±/2 and 2^Jm1 + p\ has a small effect on the predicted asymme- 
tries. The CTEQ5L [^3[ parton distribution functions are used. For this parton set, it is 
appropriate to use the one-loop running strong coupling constant with four active flavors 
and Aqcd = 0.192 GeV. Since the E691 and E687 collaborations used a beryllium target, 
the appropriate parton distribution is a linear combination of proton and neutron parton 
distribution functions: f i/Be = (4/9)/ i/p + (5/9)/ i/n . 

We use the following fragmentation function for the D meson: 

D c ^ D (z) = D(z;e)f c ^ D , (23) 

where f c ~>D is the probability for the direct fragmentation of c into D and D(z;e) is the 
Peterson fragmentation function [24]], normalized so that Jq D(z;e)dz = 1. Recent OPAL 



measurements find / c _,o»+ = 0.22 ± 0.02 p5[ and f c ->D*° = 0.22 ± 0.07 p6[. Heavy-quark 



spin symmetry implies that the direct fragmentation probabilities for D + and D° are smaller 
by a factor of 3. The probability for fragmentation to is smaller than the probability 
for D° and D + by roughly a factor of two 0. In this paper, the fragmentation probabilities 
are chosen to be 

/ c -+d*+ =3./U D+ = 0.22, (24) 
fc-^D*o = 3 f c ^ D o = 0.22 , (25) 
f^ Dt+ =3f c ^ Dt = 0.ll. (26) 

These direct fragmentation probabilities add up to about 73%, and the remaining probability 
comes from fragmentation to charmed baryons. 

Fragmentation softens the p± and xp distributions of the D mesons relative to the c 
quarks produced in the hard scattering process. The measured value of the Peterson pa- 
rameter e depends on the order to which the perturbative factors have been calculated. The 
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measured values of e corresponding to leading order perturbative calculations range from 



0.04 to 0.13 depending on the D hadron |22]. It is worth noting that the measured values 



of e for D mesons tend to be larger than for D* mesons. Larger e corresponds to a softer D 
meson distribution. The softening of the D meson fragmentation function is due primarily 
to the feeddown contributions from D* decay. For simplicity, we neglect the softening of 
the D meson fragmentation function and use a common value e = 0.06 for all D and D* 
mesons. Varying e between 0.04 and 0.13 has a minimal effect on the predictions for the 
asymmetries. 

Measurements of the charge asymmetry for a specific D meson are usually expressed 
either in terms of the anticharm/charm production ratio defined by 

R[D] = (2?) 

or in terms of the asymmetry variable 

a[D) = ^j-^ . (28 ) 

The E687 and E691 measurements of R[D] for different species of D mesons are collected 
in Table |[ The E687 measurements of a[.D + ] as a function of E^, p\ and xp are shown 
in Figs. 0, ^ and |J respectively. The E687 results are not corrected for the acceptance 
from a cut on the number of charged tracks at the photon interaction vertex, because the 
E687 hadronization model introduced a strong correlation between the acceptance and the 
asymmetries. This introduces an unknown systematic error in the comparison with theory. 
For both E687 and E691, we have combined the results from the 2 different decay channels 
of D* + into single numbers. The E687 result for the 2 different decay channels for D° can be 
combined into the single number R = 1.035±0.021. Unfortunately, the D° sample from E687 
was subjected to a "no D* tag" requirement that reduced the feeddown from D* + — > D° 
to about 20% of the total without reducing the feeddown from D*°. We have attempted to 
undo the effects of that requirement, which only complicates the comparison with theory. A 
naive prediction for the fraction of D 0, s from D* + feeddown can be obtained by neglecting 
heavy-quark recombination and using isospin and heavy-quark spin symmetry in Eq. (plj) : 
3-0.677/ (1+3-0.677+3) = 34%. A feeddown of 20% can be obtained by reducing the fraction 
multiplying a[D* + ] in Eq. (0) to 0.333: 3-0.333/(1 + 3-0.333 + 3) = 20%. The inclusive D° 
events should therefore consist of approximately (1 + 3 • 0.333 + 3)/(l + 3 • 0.677 + 3) = 83% 
D° with no D* tag events and 17% D* + events. Since the observed values of a are small, the 
value of a[_D°] is simply given by the appropriately weighted values of a[-D°,no D* tag] and 
a[-D* + ], up to corrections of order a 2 . We therefore take the E687 result for the asymmetry 
variable a for D° to be 

a[D°] =0.83a[L> ,noL>* tag] + 0.17 a[D* + ] . (29) 

We obtain a[D°] = —0.023 ± 0.009, which corresponds to the value of R listed in Table [IJ 

The E687 experiment has a broad-band photon beam with average energy (Ey) = 
200 GeV. To calculate the p]_ and xp distributions, theoretical cross sections were convoluted 
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with the photon beam spectrum of the E687 collaboration.]] Our results are nearly identical 
if instead of using the real spectra, a monochromatic 200 GeV beam is used instead. For 
the E 1 distribution, the data correspond to events collected in energy bins centered around 
40, 120, 200 and 280 GeV. For simplicity, we calculated each data point with fixed photon 
energies corresponding to the center of each bin. Feynman xp is defined by 

x F = * * m ^ Sinh( ^ , (30) 

\Fz | max 

where = p 2 ± + m 2 , y is the rapidity of the D in the photon- nucleon center-of-mass 
frame and s^n is the photon-nucleon center-of-mass energy squared. The positive z-axis is 
defined by the incident photon direction. We use the approximate expression for xf in our 
calculations. The experimental results were all subject to the cut xp > 0. 

To fix the parameters p sm and p s f, we did a least-squares fit to the E687 data consisting 
of the ratios R[D] for the four mesons in Table H, the p 2 ± distribution for a[.D + ] consisting 
of the five data points in Fig. ^|, and the xf distribution for a[-D + ] consisting of the first 
four data points in Fig. f|. Data from the last xf bin was not included in the fit because the 
theory is expected to be less reliable near xp = 1. In this corner of phase space, one expects 
increased uncertainty due to higher order corrections, nonperturbative fragmentation effects, 
the difference between the charm quark mass and D meson mass, etc. The results of a two 
parameter fit gives a value of p s f that is small and negative. Since p s f is proportional to a 
fragmentation probability, this negative value is unphysical. Therefore, we set p s f = 0, and 
the fit then yielded p sm = 0.15. 

Theoretical predictions for the kinematic dependence of a[.D + ] are compared to the E687 
results in Figs. 0, |3| and f|. Heavy-quark recombination describes the observed asymmetries 
well. In particular, the E 1 dependence and the p\ distribution are fit very well. The 
agreement with the xp distribution is not as good, but the largest discrepancy is about 2.5 
a in the last xp bin, where theoretical errors are greatest. Predictions for the R[D] ratios 
using the parameters p sm = 0.15 and p s f = are compared to the E687 and E691 data in 
Table |I[ Most predictions for the R[D] ratios are within 1 a of the data. However, the 
predictions for the R[D* + ] ratios of both experiments are below the data by about 2.6 a. A 
better fit to the R[D] ratios could be obtained if the larger value p sm = 0.3 is chosen, but 
at the expense of a worse fit to the kinematic distributions for a[.D + ]. 

The heavy-quark recombination mechanism reproduces some interesting qualitative fea- 
tures of the data. One feature is that the asymmetries decrease as E 1 increases. There are 
several effects that all tend to decrease the asymmetry. First, the ratio of quarks to gluons 
at small x decreases with energy. Furthermore, the ratio of q to q at small x approaches 
1, further diluting the asymmetry. Finally, the parton cross sections for heavy-quark re- 
combination fall faster than photon-gluon fusion as the parton center-of-mass energy grows. 
A second qualitative feature is that R[D + ], R[D* + ] > R[D ]. This is largely because the 
recombination cross section for cd is greater than that for cu. Because the beryllium nu- 
cleus is very close to an isosinglet target, the difference in these cross sections is not due 



1 A convenient parametrization of this spectrum can be found in Section 5 of 
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to the parton distribution functions, but rather to the parton cross sections in Eq. (|6]) and 
Eq. (^). The differences appear because the coupling to the photon does not respect isospin 
symmetry: k = 1 for cu and K = — | for cd. 

In Fig. |5|, the predicted inclusive cross sections, da/dxp, for D + and D~ mesons are 
shown. Also shown is the contribution to the cross section from the leading order photon- 
gluon fusion process. Again we have used the E687 photon beam spectrum with (E 7 ) = 
200 GeV. For p sm = 0.15 and p B f = 0, heavy-quark recombination contributes about 25% of 
the total cross section for both the D + and D~ . The cross sections are small for xp < 0, 
however, in this region the D~ cross section is much greater than the D + cross section. 
Thus, heavy-quark recombination predicts that |a[-D + ]| is largest in the region of phase 
space excluded by the experimental cut xf > 0. This is because recombination is most 
important when the (cq) emerges in the forward direction of the light quark in the nucleon, 
as shown in Eq. ([TT|) . 

The calculations of this paper neglected NLO corrections to the leading order photon- 
gluon fusion diagrams. If the NLO corrections are incorporated by including a K factor, 
the extracted value of p sm must be multiplied by the same K factor to obtain the same 
asymmetry. As stated earlier, there is at least a 30% uncertainty in this parameter due to 
1/rriQ, £77(3) breaking and 1/N C corrections, and a similar size uncertainty in p s f. 



IV. OTHER APPLICATIONS 



The heavy-quark recombination mechanism also can be used to calculate the asymme- 
tries of charmed baryons such as the A+. Although the cq recombination mechanism will 
contribute at some level, the most important contribution is expected to come from cq re- 
combination. We expect the asymmetry to be dominated by the analog of process a) in 
Eq. (|T8|), which would lead to a positive asymmetry a[A+]. In the measurements from E687 
and E691 in Table I, the central values correspond to a positive asymmetry, but the errors 
are large enough that the asymmetry is not statistically significant. The FOCUS experiment 
at Fermilab will measure the asymmetries for A+ with much higher statistics and will also 
measure the asymmetries for other charmed baryons, such as S°, E+ + and It will be 

interesting to compare these measurements with the predictions of our heavy-quark recom- 
bination mechanism. Quantitative predictions will involve nonperturbative factors of the 
form p[{cq) n — > A+]. Once these are fit to reproduce the overall asymmetries, the kinematic 
distributions of the asymmetries will be predicted. 

The heavy-quark recombination cross sections calculated in this paper predict that ct[/J+] 
should have the opposite sign as a[/J ]. However, it is important to point out that 

the — D~ asymmetry is generated by the fragmentation of charm quarks from process 
b) in Eq. (|l^) in which the charm antiquark is involved in the recombination process. In cq 
recombination, there is a process analogous to process b) which can dilute the asymmetry 
generated by cq recombination. To provide reliable predictions for «[-Dj~] asymmetries, both 
cq recombination into mesons and cq recombination into baryons need to be included. 

Recent photoproduction measurements by the ZEUS detector at HERA seem to indicate 
that NLO calculations underestimate the D* production rate in the forward proton direction 
1 27 1, although this effect has not been confirmed by the HI detector |p8[| . If there is an excess 
of D* in the forward direction, our heavy-quark recombination mechanism may be able to 
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explain it. Heavy-quark recombination effects should be most important in the forward 
proton direction, because the ratio of the partonic cross sections for cq recombination and 
photon-gluon fusion is highest in this direction. The forward nucleon direction is excluded 
in fixed target photoproduction measurements by the xp > cut, but it is experimentally 
accessible at the HERA ep collider. For quantitative calculations of the charm asymmetries 
at HERA energies, it will be necessary to take into account the contribution of resolved 
photons. 

A large charm asymmetry has been observed in fixed-target hadroproduction experi- 
ments, where it is called the "leading particle effect" [|29[| : charmed hadrons sharing a valence 
quark in common with the beam hadron are produced more copiously than other charmed 
hadrons in the forward region. Our heavy-quark recombination mechanism can be used to 
predict these asymmetries. The relevant parton processes are the same as those considered 
in [JTTJ] for B meson production at the Tevatron. Because fixed-target hadroproduction in- 
volves lower energies, heavy-quark recombination is expected to play a more important role 
in these experiments. It would be interesting to see if heavy-quark recombination could 
provide a quantitative description of the leading particle effect. 

We thank P. Garbincius, R. Gardner and J. Wiss for providing useful information about 
the E687 experiment. This work was supported by the Department of Energy under Grant 
No. DE-FG02-91-ER4069 and by the National Science Foundation under Grant No. PHY- 
9800964. 
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TABLE I. Anticharm/charm production ratio R = Njj/Nd: experimental results from E687 
and E691, and theory predictions with p sm = 0.15, p s f = 0. 



Particle 


E687 § 


Theory 


E691 @ 


Theory 




(Ey) = 200 GeV 




(E y ) = 145 GeV 




D+ 


1.08 ±0.02 


1.06 


1.04 ± 0.03 


1.07 


D° 


1.047 ±0.018 


1.03 


1.08 ±0.03 


1.03 


D*+ 


1.112 ±0.024 


1.05 


1.190 ± 0.049 


1.06 


Df 


0.95 ±0.10 


0.91 


0.92 ±0.14 


0.88 


A+ 


0.93 ±0.14 




0.79 ±0.17 




D*° 




1.02 




1.02 


D* + 

^ s 




0.91 




0.88 
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FIG. 2. The dependence of a[-D + ]. The data points are the measurements of E687 and 
the histogram is our prediction with p sm = 0.15, p s { = 0. 
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FIG. 3. The p\ distribution of a[L> + ]. The data points are the measurements of E687 and 
the histogram is our prediction with p sm = 0.15, p s f = 0. 
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FIG. 4. The xj? distribution of The data points are the measurements of E687 and 

the histogram is our prediction with p sm = 0.15, p s f = 0. 
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Abstract 



The asymmetries between charm and anticharm mesons observed in fixed- 
target photoproduction experiments are an order of magnitude larger than 
the asymmetries predicted by conventional perturbative QCD. We show that 
these charm meson asymmetries can be explained by a heavy-quark recom- 
bination mechanism for heavy meson production. In this process, a charm 
quark combines with a light antiquark from the hard-scattering process and 
they subsequently hadronize into a state including the charm meson. This 
recombination mechanism can be calculated within perturbative QCD up to 
some nonperturbative constants. After using symmetries of QCD to reduce 
the number of free parameters to two, we obtain a good fit to all the data on 
the asymmetries for charmed mesons from the E687 and E691 experiments. 



Typeset using REVTeX 
*Address after Jan. 1, 2002: Department of Physics, Duke University, Durham, NC 27708. 



High energy charm photoproduction is an excellent testing ground for our understand- 
ing of strong interactions Q. The charm quark mass provides a large scale justifying the 
use of perturbative QCD. Theoretical uncertainties in photoproduction are smaller than in 
hadroproduction because only one hadron is present in the initial state. Charm events also 
constitute a higher fraction of the total hadronic cross section than in hadroproduction, 
which facilitates experimental studies. 

High energy charm photoproduction experiments [@-|J] have observed an interesting pro- 
duction asymmetry between charm and anticharm mesons. The E691 and E687 Q 
experiments at Fermilab measured the cross sections for various D and D mesons from a 
photon beam on a beryllium target. The experiments imposed a cut on the Feynman xp 
variable, xp > 0, which corresponds to the forward direction of the photon beam in the 
7 A?" center-of-mass frame. E691 and E687 found statistically significant asymmetries in 
D+-D-, D°- D and D*+ - D*~. In each case the cross section for the anticharm meson 
is greater, though the size of the asymmetry is different for each species of meson. Produc- 
tion asymmetries for and A+ are also measured, but with much larger errors so they 
are consistent with zero. For D + and D~ mesons, the E687 experiment has also measured 
the dependence of the charm asymmetry on the photon energy, E 11 the meson transverse 
momentum, p±, and Feynman xp. 

The observed charm asymmetries in photoproduction are an order of magnitude larger 
than the asymmetries between charm and anticharm quarks predicted by perturbative QCD. 
Thus they cannot be understood within the fragmentation mechanism for heavy meson pro- 
duction which is based on the factorization theorem for inclusive single particle production 
0. In this mechanism, the c and c are produced in a hard-scattering process and then 
fragment independently into charmed hadrons. At leading order in a s , charm is produced 
by photon-gluon fusion |J which produces c and c symmetrically. The next-to-leading order 
(NLO) correction increases the normalization of the total cross section by about 50% 
at the energies relevant to the E687 and E691 experiments. The shapes of the xp and p± dis- 
tributions are similar to those at leading order, but there is a tiny charm asymmetry in these 
kinematic distributions, arising from the processes 7 + q(q) — > c + c + q(q). For E 1 = 200 
GeV, the NLO calculations predict the ratio of the c and c cross sections in the forward 



region to be about R = 1.006 jlOfl . The asymmetry, (1 — R)/(l + R) = —0.003, is about an 
order of magnitude smaller than the asymmetries for D mesons measured by the E687 and 
E691 experiments. Furthermore, the NLO calculations cannot explain the differences in the 
asymmetries of D mesons with different light quark flavors. It is commonly assumed that a 
nonperturbative hadronization mechanism involving the remnant of the nucleon or photon 
after the hard scattering is responsible for the charm asymmetries. 

In this paper, we argue that the charm asymmetries can be explained naturally and 
economically by the heavy-quark recombination mechanism for heavy meson production in- 
troduced in In the cq recombination mechanism, a light q from the target nucleon 



participates in the hard-scattering process that produces the c and c. The light q emerges 
from the hard-scattering with small momentum in the c rest frame, and the c and q subse- 
quently hadronize into a final state including the D meson. There is also a cq recombination 
process in which a light q from the target nucleon recombines with a c to hadronize into a 
D meson. Because the parton distribution functions of u and d in the nucleon are greater 
than the parton distribution functions for u and d, the recombination cross sections for cu 



and cd are greater than for cu and cd. This gives rise to the asymmetries between charm 
and anticharm mesons. Since the photon participates in the hard scattering and can couple 
to the light quark, the recombination cross sections for processes involving u quarks differ 
from those involving d and s quarks because of their different electric charges. Thus, the 
heavy-quark recombination mechanism can account for the differences in the asymmetries 
of D mesons with different light quark flavor quantum numbers. Below we will see that it 
can also adequately describe the kinematic dependence of the asymmetries on E^, p± and 

Xp. 

Before discussing the heavy-quark recombination mechanism in detail, we briefly review 
the models which have been proposed in the literature to explain the D meson asymmetries. 
Though these models differ in details, all generate the asymmetries by a nonperturbative 
hadronization process that involves the remnant of the nucleon or photon after the hard 
scattering. The predictions of these models depend on unknown and ad-hoc functions that 
specify the distribution of partons in the remnant. 

The model for the asymmetry that has been most thoroughly developed involves the 
fragmentation of strings connecting the c and c to the nucleon remnant [|l2j. After the c 
and c are produced via photon-gluon fusion, the remnant of the target nucleon is left in a 
color-octet state. The remnant is modelled as a color-triplet "bachelor quark" and a color- 
antitriplet diquark with momentum fractions x and 1 — x specified by a distribution function, 
F(x). The subsequent hadronization is calculated using the Lund string fragmentation 
model Color-field strings are introduced between the diquark and the c and between 
the bachelor quark and the c. The prediction for the asymmetry depends sensitively on the 
choice of F(x) as well as the relative probability for each string to break into a charmed 
meson or a charmed baryon. This model reproduces the qualitative shapes of the kinematic 
distributions for the D + — D~ asymmetry measured by the E687 collaboration ||]. However, 
the magnitude of the asymmetries is sensitive to the unknown function F(x). For a "soft" 
distribution that diverges as 1/x at small x (the default option in the PYTHIA 5.6 event 
generator program), the model tends to overestimate the asymmetry in all kinematic regions. 
However, the data can be fit with a "hard" momentum distribution proportional to 1 — x [Q. 

Models for the asymmetry that involve the remnant of the hadronic component of the 
photon have also been proposed. In one such model a c is produced by the hard 
scattering process qq — > cc involving a q from the photon. The c then recombines with a q 
from the photon remnant to form a D meson. The predictions depend on the distribution 
function for the q in the photon remnant and on a recombination function that could be 
calculated in principle using the Lund string fragmentation model. Another model is based 



on the possibility of intrinsic charm JT3J in the photon. In this model [IS], a q in the resolved 



photon undergoes a hard scattering and the D forms from a c and q in the photon remnant. 
The prediction depends on the distribution function for the momentum fractions of the c and 
q in the photon remnant and a recombination function that depends on several variables. 
Both of the models involve arbitrary functions that can be tuned to fit the data. 

In contrast with these models, the heavy-quark recombination mechanism [] L| generates 
the asymmetries within a hard-scattering process. The asymmetries are calculable in per- 
turbative QCD up to some nonperturbative parameters related to the probability for the c 
and q produced in the short- distance process to hadronize into a D meson. Because the c 
and q do not necessarily have the same color, angular momentum or light quark flavor quan- 
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turn numbers as the final state D meson, several nonperturbative parameters appear in the 
calculation. However, it is possible to use heavy quark symmetry, 577(3) flavor symmetry 
and large N c arguments to reduce the number of parameters to two. The asymmetries can 
be adequately described by fitting these two parameters. Because of this, the heavy-quark 
recombination mechanism gives a very economical description of the observed asymmetries 
and does not rely on model-dependent assumptions about nonperturbative hadronization. 

The factorization theorems of QCD || show that the leading contribution to the D 
meson photoproduction cross section at large transverse momentum can be written as: 

do^i + N -> D + X] = fi/N ® da[y + i -f c + c + X] <g> D-^ , (1) 

i 

where fam is the parton distribution function for the parton j in a nucleon and D-_^ is 
the nonperturbative fragmentation function for c to hadronize into D. At leading order 
in a s , the only parton that contributes is the gluon, so photon-gluon fusion dominates the 
cross section. At large p±, Eq. ([]]) receives corrections suppressed by powers of Aq C d/p± or 
m c/p±- The heavy-quark recombination mechanism introduced in [[Uj is an 0(AqcD m c/ Pj_) 
power correction to Eq. ([I]). The eg recombination contribution to D production can be 
expressed as 

da[y + N^D + X] = f q/N ®J2 da h + ( l^ (c<?) n + c] p[(cq) n -> D] . (2) 

n 

Making the replacement c <-> c, q — > g and U — > JJ in Eq. (Q) gives the eg recombination 
contribution to D production. The notation (cg) n in Eq. @ indicates that the g in the 
final state has momentum of 0(Aq C d) in the c rest frame, and that the c and the g are in 
a state with definite color and angular momentum quantum numbers specified by n. The 
factor p[(cq) n — >• D] is a nonperturbative constant proportional to the probability for (cg) n 
to evolve into a final state that includes the D. The recombination process in Eq. (0) is an 
0{AQCD m c/v\) power correction to Eq. ([!]), because the parton cross section in Eq. @ is 
suppressed by 0{m 2 c /p\) relative to the parton cross section in Eq. (|I|) and p[(cq) n — > D] is 
expected to be 0(Aqcd/^c) from heavy quark effective theory scaling arguments. Because 
the momentum of the light quark is O(Aqcd), the production of (cq) n in states with L / 
is suppressed by powers of Aq C d/p± or Aqcd/^c relative to S-waves. Therefore, it is only 
necessary to compute production of (eg)™ with 1 S'o or 3 Si angular momentum quantum 
numbers. The (cq) n can be in either a color-singlet or a color-octet state. 

Berezhnoy, Kiselev, and Likhoded ]TB|-p^|] have developed a model for heavy meson 
production that incorporates some aspects of our heavy-quark recombination mechanism. 
Their model is based upon an extrapolation of the NRQCD factorization formalism for heavy 
quarkonium production ||20|| . A D or D* meson is modelled as a nonrelativistic quarkonium 
composed of a c and q with masses m c and m q . They take m q = 300 MeV, consistent with an 
interpretation as a constituent quark mass. Their expression for the cross section is obtained 
by keeping only specific terms in the NRQCD factorization formula. These terms involve 
S-wave matrix elements, which are distinguished by the fact that their contributions at large 
p± are suppressed only by m q m c /p\. Of the four independent S-wave matrix elements, they 
keep only two: the Color-singlet Model matrix element (Oi( 3 Si)) D * = 3((9i( 1 S'o)) i:) and the 
color-octet matrix element {Os( 3 Si)) D * = 3(Os(}So)) D . Thus their model involves only 3 
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FIG. 1. Feynman diagrams for the recombination process 7 + q — > (cq) n + c. Single lines 
represent a light quark and double lines represents a charm quark. The shaded blob represents 
(cq) n . 

parameters: m g , (Oi^Sq)) and ((9 8 ( 1 5'o)) D . In their model, the terms in the cross sections 
that diverge like l/m q as m q — > have the same form as our heavy-quark recombination 
cross sections if we make the following identification between the nonperturbative factors: 
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Our approach provides a model-independent framework for calculating those contributions. 



In the model of Ref. fl§ |, contributions corresponding to the NRQCD matrix elements 
(Oi( 3 Si)) D and (C 8 ( 3 5'i)) r) , which require a spin-flip between the partons and the meson, 
are omitted. In the heavy quarkonium case, these matrix elements are suppressed by v 7 and 
v 3 with respect to (Oi^So)) , while (Os^So)) is suppressed by t> 4 . But when the antiquark 
q becomes light, v approaches 1 and there is no parametric suppression that can be used to 
justify omitting them. Furthermore, the NRQCD operators Ci( 1 S , ), Ci( 3 «Si), 8 ( 1 5 , o ), and 
0$( 3 Si) mix under renormalization. This implies that at some order in perturbation theory, 
inconsistencies will arise if any of these operators are omitted. It is possible that some of 
them have matrix elements that are numerically small, but they must all contribute at some 
level. 



The most important difference between our approach and the model of ref. |19| is that 
we use our nonperturbative probability factors p only for very specific terms in the cross 
sections: those that diverge as l/x q , whre x q is the light-cone momentum fraction of the q 



relative to the heavy quark. In contrast, the model of Ref. |19[ is applied in all kinematical 
regions. Their parameters m q , (Oi( 1 5'o)) D , and (Os^So)) are used to calculate not only 
the heavy-quark recombination contribution to the cross section, but all other contributions 
as well, such as the fragmentation function D c ^d(z). In contrast, we take the fragmentation 
function to be an nonperturbative function that is independent of the probability factors 



p for heavy-quark recombination. We use fragmentation functions which have been been 
extracted from D meson production in e + e~ collisions. With only 3 parameters, the model 
of Ref. [B| is very constrained, but the predictive power comes at a sacrifice in rigor. 



The four Feynman diagrams that contribute to the recombination process 7+g — > (cq) + c 
at lowest order in a s are shown in Fig. [I]. Since the momentum of the light quark in the final 
state, p' , is O(Aqcd) in the c rest frame, its transverse momentum can be neglected and 
its light-cone momentum fraction, x q , is 0(AQCD/ m c)- in the limit x q — > 0, the diagrams in 
Figs. 0a), b) and d) are 0(x~ 1 ) while the diagram in Fig. |l|c) is 0(x° q ). Therefore the diagram 
in Fig. |T|c) gives an 0(Agc£)/m c ) suppressed contribution to the heavy-quark recombination 
rate and can be neglected. The calculation is very similar to the calculation of the process 
g + q — > (bq) + b in Ref. Interested readers can refer to [JTTJ] for a detailed description of 



the calculation. The color-singlet 1 S'o contribution to the cross section is obtained by first 
making the following substitution in the parton amplitude: 

Vi(Pc)uj(p q ) — ► x q j^-m c f + (fc - m c )7 5 . (4) 

We then set p' q = x q Pc in the rest of the amplitude and take the limit x q — * 0. The 
nonperturbative parameter / + is related to the amplitude for the (cq) produced in the short- 
distance process to bind into a D meson. /+ is proportional to the D meson decay constant, 
fo, times a moment of the D meson light-cone wave function, and therefore is expected to 
scale as m" 1 / 2 in the heavy quark limit. However, / + does not account for the possibility 
that the (cq) can hadronize into states that include other light hadrons whose momenta are 
soft in the D meson rest frame. In order to take these final states into account, j\ is replaced 
with the parameter pi[cg( 1 S'o) — > D] after squaring the matrix element. Pi[cg( 1 S'o) — > D] is 
proportional to the inclusive probability for the color-singlet, spin-singlet (cq) to evolve into 
any state containing the D hadron. The calculation of the color-singlet 3 Si cross section is 
identical except that 75 in Eq. (Q) is replaced by /, where e M is a polarization vector, and 
f\ is replaced by pi[cq( 3 Si) — ► D}. 

The D production cross section also receives contributions from processes in which the 
(cq) is produced in a color-octet configuration in the short-distance process. The color-octet 
1 5'o cross section is obtained by making the substitution 

Vi(Pc)uj(p q ) — ► XqJ^-Tfj m c fl(0c - m c )7 5 (5) 

in the parton amplitude, replacing p' — > x q Pc and then taking the limit x q — > 0. After 
squaring the matrix element and summing over colors in the final state, (/+) 2 is replaced 
with p 8 [cg( 1 5'o) — > D]. The normalization in Eq. @ is chosen so that ps/pi is the ratio of 
the probabilities for a color-octet (cq) and a color-singlet (cq) to hadronize into a D hadron. 
Finally the color-octet 3 Si contribution is obtained by replacing 75 in Eq. @ by and 
replacing (/|) 2 by p 8 [cq( 3 Si) -> D]. 

The final results for the parton cross sections for color-singlet cq recombination in 7 + q 
collisions are 
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where k = e q /e c is the ratio of the electric charge of the light quark q to that of the charm 
quark, S = (p g + p 7 ) 2 , T = (p 7 — p c ) 2 — m 2 and {7 = (p 7 — Pc) 2 — r^ 2 . These variables are 
defined so that S 1 + T + [/ = 0. If ^ is defined as the angle between the incoming q and the 
outgoing (cq) in the parton center-of-mass frame, then T = —S/2[l — (1 — Am 2 J Sf 2 cos 9]. 
The contributions to the cross section for D meson production are obtained by multiplying 
the parton cross sections in Eqs. flf]) and (0) by pi[cq(}S ) — ► Z5] and p 1 [cg( 3 5'i) — > D], 
respectively. Note that the recombination process is calculated in the heavy quark limit, so 
we do not distinguish between m c and the D meson mass. 

For the photoproduction diagrams in Fig. 0, the color-octet amplitude is identical to 
the color-singlet amplitude up to an overall color factor. The color-octet cross sections can 
be obtained from the color-singlet cross sections in Eqs. (|) and (0) by replacing p\ with 
Ps/8. Therefore, heavy-quark recombination contributions to photoproduction cross sections 
depend on the following linear combination of color-singlet and color-octet parameters: 



PcS = Pi + o PS 



(8) 



It is interesting to compare the relative sizes of the heavy-quark recombination and 
photon-gluon fusion cross sections in different regions of phase space. Setting 9 = n/2 and 
expanding in powers of m 2 /S, one finds that the ratio of the parton cross sections is 
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(10) 



For 9 = n/2, S = 4:(p 2 ± + m 2 ), so the parton cross section for heavy-quark recombination 
is suppressed by ml/p 2 ± at large The recombination cross sections are therefore sup- 
pressed by Aqcd'>tT'c/p± relative to the photon-gluon fusion cross section, in accord with the 
factorization theorem. The differential cross section is finite as 9 — > and 9 — > it because 
the heavy quark mass acts as an infrared cutoff. Setting 9 = and expanding in rn^/S, one 
finds that the ratio of the two parton cross sections is 



da[y + q — > (cq) + c] 



da[y + g — > c + c] 



256 7T 



for cqCS ),cq( 3 Si). 



(11) 



Thus, the heavy-quark recombination cross sections have no kinematic suppression factor 
when the (cq) emerges from the parton collision in the same direction as the incident q. On 
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the other hand, the recombination cross sections are suppressed by factors of m 2 c /S when 
the (cq) emerges in the same direction as the incident photon: 



da[y + q — > (cq) + c] 



da[y + g — > c + c] 



25 p« s || for cqCSo), (12) 

\ ± ^a s ^ for cq( 3 Sx) . (13) 



e=7r 



The kinematic suppression is especially strong in the 1 Sq case. 

To analyze the photoproduction data of E687 and E691, it is necessary to compute the 
cross section for various species of D mesons. It is important to realize that the color, 
angular momentum and light quark flavor quantum numbers of the (cq) n are not necessarily 
the same as the D meson in the final state [O. The transition from the (cq) n to the final 



state that includes the D meson is governed by nonperturbative physics that can change 
these quantum numbers. For instance it is possible for the (cq) to be produced in a color- 
octet state in the short-distance process and emit soft gluons during the nonperturbative 
transition to emerge as a color-singlet D meson. It is also possible for the spin of the light 
quark to be flipped during this transition, resulting in a D meson with different angular 
momentum. Finally, a light q'q' pair can be created in the transition, followed by the c 
binding with the q' to form a D meson with different flavor quantum numbers than the 
original (cq) n . Thus for a specific D meson, several nonperturbative parameters enter the 
calculation. For instance, for the D° meson, the parameters are p e fi[cu( l So) — > D°] and 
p e s[cu( 3 Si) — > D°], together with those obtained by replacing u with d or s. 

Symmetries of the strong interaction allow us to greatly reduce the number of free pa- 
rameters. Heavy quark spin symmetry |^T| implies that the spin of the c quark should 



remain unchanged throughout the nonperturbative transtion. This means that transitions 
which can be related by flipping the spin of the heavy quark in both the initial and final 
state should have the same rate. For example, 

PoslcdCSo) -> D°] = p eff [cd( 3 5 1 ) -> D*°] . (14) 

One can also use SU(3) light quark flavor symmetry to relate nonperturbative transitions 
with different light quark flavors. An example of an SU(3) relation is 

p cS [cuCS ) -> D°] = PeS [cdCSo) - D + ] = p cS [csCS ) -> £>+] . (15) 

Finally, one can argue that in the large iV c limit of QCD, light quark-antiquark pair produc- 
tion is suppressed. This suppression in low energy QCD is the basis for the phenomenolog- 
ically successful Zweig's rule. Suppression of light quark-antiquark pair production in the 
long-distance transition then forces the (cq) n to have the same flavor quantum numbers as 
the final state D meson. 

Using the symmetries discussed above together with the large N c argument, one is left 
with two parameters. One is 

Psm = Pcs[cuCS ) -> £°] = pcslcdCSo) -> D+] = p eS [cs( l S ) -> D+] 

= PeS [(mCSi) - = peffK 3 ^) -> d* + ] = PeS [c8( 3 Si) - d* s + ] . (16) 
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The subscript sm stands for spin-matched, since in all these transitions, the spin of the cq 
is the same as the spin of the D meson. The other parameter describes transitions in which 
the light quark has its spin flipped: 

Psf = Pos[cu( 3 Si) -> D°] = pefrM 3 ^) -> D+] = p eff [cs( 3 £i) -> D+] 

= PcS [cuCS ) - 7J*°] = Peff K 1 ^) - D*+] = prticsCSo) - . (17) 

While the use of heavy quark and £77(3) symmetry greatly reduces the number of parame- 
ters, it should be kept in mind that these relations are only expected to hold at the 30% level. 
It is also important to note that the flavor-changing transitions are not strongly suppressed 
because N c = 3. The relevant nonperturbative parameters for these transitions are expected 
to be smaller than p sm and p s f but nonvanishing. These are not included in the calculations 
of this paper because an adequate fit to E687 and E691 data can be obtained using only 
the two parameters p sm and p s f. However, it may be necessary to include flavor-changing 
transitions in other applications, or if the accuracy of current data on charm asymmetries 
significantly improves. 

It is important to note that when a D is produced by a cq recombination process, 
the associated c quark can produce a D meson via the usual fragmentation mechanism. 
Therefore, the recombination process can contribute to the direct cross section for D mesons 
in two ways: 

a ) £ daYl + q^ (cq) n + c] p n [(cq) n - D] , (18) 

n 

b) d °h + <l ( s 9) n + c] Pn ® D C ^ D . (19) 

q, n 

In a), the (cq) n recombines into the D meson. In b), the (cq) n recombines into a D meson, 
and the D comes from the fragmentation of the recoiling c. For direct D production, Eq. (|I8"D 
and Eq. (|19"D are replaced with their charge conjugates. Process a) produces more D~ and 
D° mesons than D + and D°, simply because the nucleon contains more u and d than u 
and d. In process b), D's are produced when there is a q in the initial state, while JJ's are 
produced when there is a q in the initial state. The excess of u and d over u and d in the 
nucleon then leads to more D than D. It turns out that in the case of D + — D~ and D° — D 
mesons, process b) dilutes the asymmetry produced in a). For D s mesons, process a) cannot 
generate an asymmetry because the s and the s content of the nucleon are identical. Process 
b) generates an asymmetry for Df — D~ which has the opposite sign as the D + — D~ and 
D° — D asymmetries. 

Finally, when computing the inclusive cross section for D and D* mesons, one must 
consider the possibility of feeddown from other D mesons that are produced via photon-gluon 
fusion or heavy-quark recombination and then decay to D or D*. The feeddown from excited 
D mesons, e.g. is not included but is expected to be small. Therefore, the inclusive 
D* cross section is the sum of the photon-gluon fusion and heavy-quark recombination 
contributions, which will be referred to as the direct cross section. The inclusive D meson 
cross section also includes feeddown from D* mesons. Denoting the direct cross sections by 
o"dir an d the inclusive cross sections by <j inC) one finds 

a inc [D + ] = a dil [D + ] + 0.323 a dir [/J*+] , (20) 
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<Xinc[£°] = a dir [D°] + 0.677 a dil [D* + ] + a dir [D*°] , (21) 
a inc [Df] = a dil [Df] + a dir [D* s + ] . (22) 



The coefficients of <7di r [-D* + ] in Eq. (20) and Eq. (|21|) are the branching fractions of D* + into 
D + and D° | 22| . The cross sections for the charge conjugate states have similar expressions. 

In the calculations presented in this paper, the contribution to charm photoproduction 
from the hadronic component of the photon is neglected because it constitutes less than 
5% of the cross section for the fixed-target experiments in the energy range 50 GeV < 
E 1 < 400 GeV 0. The charm quark mass is taken to be 1.5 GeV and both the factoriza- 



tion and renormalization scales are set equal to Jm? c + p\. We checked that varying this 



scale between ym;? + p\j1 and 2ym^ + p\ has a small effect on the predicted asymme- 
tries. The CTEQ5L [^3[ parton distribution functions are used. For this parton set, it is 
appropriate to use the one-loop running strong coupling constant with four active flavors 
and Aqcd = 0.192 GeV. Since the E691 and E687 collaborations used a beryllium target, 
the appropriate parton distribution is a linear combination of proton and neutron parton 
distribution functions: f i/Be = (4/9) f i/p + (5/9) f i/n . 

We use the following fragmentation function for the D meson: 

D c ^ D {z) = D{z;e)f c ^ D , (23) 

where f c ^D is the probability for the direct fragmentation of c into D and D(z;e) is the 
Peterson fragmentation function | 23}| , normalized so that Jq D(z;e)dz = 1. Recent OPAL 



measurements find f c ^o*+ — 0.22 ± 0.02 [25] and f c ->D*° — 0.22 ± 0.07 [26]. Heavy-quark 



spin symmetry implies that the direct fragmentation probabilities for D + and D° are smaller 
by a factor of 3. The probability for fragmentation to is smaller than the probability 
for D° and D + by roughly a factor of two |]. In this paper, the fragmentation probabilities 
are chosen to be 

f c ^ D * + =3./U D+ = 0.22, (24) 
7^.0 = 3/ M2J o = 0.22, (25) 
/ c _ d:+ =3/ c ^ s+ = 0.11. (26) 

These direct fragmentation probabilities add up to about 73%, and the remaining probability 
comes from fragmentation to charmed baryons. 

Fragmentation softens the p± and xp distributions of the D mesons relative to the c 
quarks produced in the hard scattering process. The measured values of the Peterson pa- 
rameter e range from 0.04 to 0.13 |22j depending on the D hadron. It is worth noting that 
the measured values of e for D mesons tend to be larger than for D* mesons. Larger e 
corresponds to a softer D meson distribution. The softening of the D meson fragmenta- 
tion function is due primarily to the feeddown contributions from D* decay. For simplicity, 
we neglect the softening of the D meson fragmentation function and use a common value 
e = 0.06 for all D and D* mesons. Varying e between 0.04 and 0.13 has a minimal effect on 
the predictions for the asymmetries. 

Measurements of the charge asymmetry for a specific D meson are usually expressed 
either in terms of the anticharm/charm production ratio defined by 
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R[D] = , (27) 

or in terms of the asymmetry variable 

a[D) = ^\ D ]-^\B . (28) 

The E687 and E691 measurements of R[D] for different species of D mesons are collected 
in Table |. The E687 measurements of a[.D + ] as a function of E^, p\ and xp are shown 
in Figs. 0, ^ and f|, respectively The E687 results are not corrected for the acceptance 
from a cut on the number of charged tracks at the photon interaction vertex, because the 
E687 hadronization model introduced a strong correlation between the acceptance and the 
asymmetries. This introduces an unknown systematic error in the comparison with theory. 
For both E687 and E691, we have combined the results from the 2 different decay channels 
of D* + into single numbers. The E687 result for the 2 different decay channels for D° can be 
combined into the single number R = 1.035±0.021. Unfortunately, the D° sample from E687 
was subjected to a "no D* tag" requirement that reduced the feeddown from D* + — > D° 
to about 20% of the total without reducing the feeddown from D*°. We have attempted to 
undo the effects of that requirement, which only complicates the comparison with theory. A 
naive prediction for the fraction of D 0, s from D* + feeddown can be obtained by neglecting 
heavy-quark recombination and using isospin and heavy-quark spin symmetry in Eq. ([ZT|): 
3-0.677/ (1+3-0. 677+3) = 34%. A feeddown of 20% can be obtained by reducing the fraction 
multiplying a[D* + ] in Eq. (0) to 0.333: 3-0.333/(1 + 3-0.333 + 3) = 20%. The inclusive D° 
events should therefore consist of approximately (1 + 3 • 0.333 + 3)/(l + 3 • 0.677 + 3) = 83% 
D° with no D* tag events and 17% D* + events. Since the observed values of a are small, the 
value of a[D°] is simply given by the appropriately weighted values of a[D°, no D* tag] and 
a[-D* + ], up to corrections of order a 2 . We therefore take the E687 result for the asymmetry 
variable a for D° to be 

a[D°] = 0.83 a no tag] + 0.17 a [£>*+]. (29) 

We obtain a[D°] = —0.023 ± 0.009, which corresponds to the value of R listed in Table [IJ 

The E687 experiment has a broad-band photon beam with average energy (E 7 ) = 
200 GeV. To calculate the p 2 ± and xf distributions, theoretical cross sections were convoluted 
with the photon beam spectrum of the E687 collaboration.^ Our results are nearly identical 
if instead of using the real spectra, a monochromatic 200 GeV beam is used instead. For 
the E 1 distribution, the data correspond to events collected in energy bins centered around 
40, 120, 200 and 280 GeV. For simplicity, we calculated each data point with fixed photon 
energies corresponding to the center of each bin. Feynman xf is defined by 

|Pz|max Y^yiV 



1 A convenient parametrization of this spectrum can be found in Section 5 of 
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where m\ = p]_ + m 2 ., y is the rapidity of the D in the photon- nucleon center-of-mass 
frame and s 7 n is the photon-nucleon center-of-mass energy squared. The positive z-axis is 
defined by the incident photon direction. We use the approximate expression for xp in our 
calculations. The experimental results were all subject to the cut xp > 0. 

To fix the parameters p sm and p s f, we did a least-squares fit to the E687 data consisting 
of the ratios R[D] for the four mesons in Table H, the p 2 ± distribution for a[.D + ] consisting 
of the five data points in Fig. and the xp distribution for a[.D + ] consisting of the first 
four data points in Fig. f|. Data from the last xp bin was not included in the fit because the 
theory is expected to be less reliable near xp = \. In this corner of phase space, one expects 
increased uncertainty due to higher order corrections, nonperturbative fragmentation effects, 
the difference between the charm quark mass and D meson mass, etc. The results of a two 
parameter fit gives a value of p s f that is small and negative. Since p s f is proportional to a 
fragmentation probability, this negative value is unphysical. Therefore, we set p s f = 0, and 
the fit then yielded p sm = 0.15. 

Theoretical predictions for the kinematic dependence of a[.D + ] are compared to the E687 
results in Figs. 0, |3| and |j. Heavy-quark recombination describes the observed asymmetries 
well. In particular, the _E 7 dependence and the p]_ distribution are fit very well. The 
agreement with the xf distribution is not as good, but the largest discrepancy is about 2.5 
a in the last xp bin, where theoretical errors are greatest. Predictions for the R[D] ratios 
using the parameters p sm = 0.15 and p s f = are compared to the E687 and E691 data in 
Table |I[ Most predictions for the R[D] ratios are within 1 a of the data. However, the 
predictions for the R[D* + ] ratios of both experiments are below the data by about 2.6 a. A 
better fit to the R[D] ratios could be obtained if the larger value p sm = 0.3 is chosen, but 
at the expense of a worse fit to the kinematic distributions for a[.D + ]. 

The heavy-quark recombination mechanism reproduces some interesting qualitative fea- 
tures of the data. One feature is that the asymmetries decrease as increases. There are 
several effects that all tend to decrease the asymmetry. First, the ratio of quarks to gluons 
at small x decreases with energy. Furthermore, the ratio of q to q at small x approaches 
1, further diluting the asymmetry. Finally, the parton cross sections for heavy-quark re- 
combination fall faster than photon-gluon fusion as the parton center-of-mass energy grows. 
A second qualitative feature is that R[D + ], R[D* + ] > R[D ]. This is largely because the 
recombination cross section for cd is greater than that for cu. Because the beryllium nu- 
cleus is very close to an isosinglet target, the difference in these cross sections is not due 
to the parton distribution functions, but rather to the parton cross sections in Eq. (g) and 
Eq. (0). The differences appear because the coupling to the photon does not respect isospin 
symmetry: k = 1 for cu and k = — | for cd. 

In Fig. ||, the predicted inclusive cross sections, dajdxp, for D + and D~ mesons are 
shown. Also shown is the contribution to the cross section from the leading order photon- 
gluon fusion process. Again we have used the E687 photon beam spectrum with (E^) = 
200 GeV. For p sm = 0.15 and p s f = 0, heavy-quark recombination contributes about 25% of 
the total cross section for both the D + and D~ . The cross sections are small for xp < 0, 
however, in this region the D~ cross section is much greater than the D + cross section. 
Thus, heavy-quark recombination predicts that |a[-D + ]| is largest in the region of phase 
space excluded by the experimental cut xp > 0. This is because recombination is most 
important when the (cq) emerges in the forward direction of the light quark in the nucleon, 
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as shown in Eq. (|TT|). 

The calculations of this paper neglected NLO corrections to the leading order photon- 
gluon fusion diagrams. If the NLO corrections are incorporated by including a K factor, 
the extracted value of p sm must be multiplied by the same K factor to obtain the same 
asymmetry. As stated earlier, there is at least a 30% uncertainty in this parameter due to 
1/rriQ, SU(3) breaking and 1/N C corrections, and a similar size uncertainty in p s f. 

The heavy-quark recombination mechanism also can be used to calculate the asymme- 
tries of charmed baryons such as the A+ . Although the cq recombination mechanism will 
contribute at some level, the most important contribution is expected to come from cq re- 
combination. We expect the asymmetry to be dominated by the analog of process a) in 
Eq. (|i~8|), which would lead to a positive asymmetry a[A+]. In the measurements from E687 
and E691 in Table I, the central values correspond to a positive asymmetry, but the errors 
are large enough that the asymmetry is not statistically significant. The FOCUS experiment 
at Fermilab will measure the asymmetries for A+ with much higher statistics and will also 
measure the asymmetries for other charmed baryons, such as E°, £+ + and It will be 

interesting to compare these measurements with the predictions of our heavy-quark recom- 
bination mechanism. Quantitative predictions will involve nonperturbative factors of the 
form p[{cq) n — > A+]. Once these are fit to reproduce the overall asymmetries, the kinematic 
distributions of the asymmetries will be predicted. 

The heavy-quark recombination cross sections calculated in this paper predict that a;[D+] 
should have the opposite sign as a[-D + ], a[-D ]. However, it is important to point out that 
the .D+ — D~ asymmetry is generated by the fragmentation of charm quarks from process 
b) in Eq. ([T^) in which the charm antiquark is involved in the recombination process. In cq 
recombination, there is a process analogous to process b) which can dilute the asymmetry 
generated by cq recombination. To provide reliable predictions for a[-Dj~] asymmetries, both 
cq recombination into mesons and cq recombination into baryons need to be included. 

Recent photoproduction measurements at HERA have shown that NLO calculations 
underestimate the D* production rate in the forward proton direction ||27|| . Our heavy- 
quark recombination mechanism may be able to explain this observation. Heavy-quark 
recombination effects should be most important in the forward proton direction, because the 
ratio of the partonic cross sections for cq recombination and photon-gluon fusion is highest 
in this direction. The forward nucleon direction is excluded in fixed target photoproduction 
measurements by the xf > cut, but it is experimentally accessible at the HERA ep 
collider. For quantitative calculations of the charm asymmetries at HERA energies, it will 
be necessary to take into account the contribution of resolved photons. 

A large charm asymmetry has been observed in fixed-target hadroproduction experi- 
ments, where it is called the "leading particle effect" ||28|| : charmed hadrons sharing a valence 
quark in common with the beam hadron are produced more copiously than other charmed 
hadrons in the forward region. Our heavy-quark recombination mechanism can be used to 
predict these asymmetries. The relevant parton processes are the same as those considered 
in [11] for B meson production at the Tevatron. Because fixed-target hadroproduction in- 
volves lower energies, heavy-quark recombination is expected to play a more important role 
in these experiments. It would be interesting to see if heavy-quark recombination could 
provide a quantitative description of the leading particle effect. 
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TABLE I. Anticharm/charm production ratio R = Njj/Nd: experimental results from E687 
and E691, and theory predictions with p sm = 0.15, p s f = 0. 



Particle 


E687 § 


Theory 


E691 @ 


Theory 




(Ey) = 200 GeV 




(E y ) = 145 GeV 




D+ 


1.08 ±0.02 


1.06 


1.04 ± 0.03 


1.07 


D° 


1.047 ±0.018 


1.03 


1.08 ±0.03 


1.03 


D*+ 


1.112 ±0.024 


1.05 


1.190 ± 0.049 


1.06 


Df 


0.95 ±0.10 


0.91 


0.92 ±0.14 


0.88 


A+ 


0.93 ±0.14 




0.79 ±0.17 




D*° 




1.02 




1.02 


D* + 

^ s 




0.91 




0.88 



0.1 



a 



0.0 



-0.1 



T 
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100 



200 



300 



E 7 (GeV) 



FIG. 2. The dependence of a[-D + ]. The data points are the measurements of E687 and 
the histogram is our prediction with p sm = 0.15, p s { = 0. 
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FIG. 3. The p\ distribution of a[L> + ]. The data points are the measurements of E687 and 
the histogram is our prediction with p sm = 0.15, p s f = 0. 
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FIG. 4. The xj? distribution of The data points are the measurements of E687 and 

the histogram is our prediction with p sm = 0.15, p s f = 0. 
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